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in 
o 
o 



Chad R. Galley and B. L. Hu 
Department of Physics, University of Maryland, College Park, Maryland, 20742 

(Dated: May 16, 2005) 

We give a quantum field theoretical derivation of the scalar Abraham- Lorentz-Dirac (ALD) equa- 
tion and the self-force for a scalar charged particle interacting with a quantum scalar field in curved 
spacetime. We regularize the causal Green's function using a quasi-local expansion in the spirit of 
effective field theory and obtain a regular expression for the self-force. The scalar ALD equation 
obtained in this way for the classical motion of the particle checks with the equation obtained by 
Quinn earlier p|. We further derive a scalar ALD-Langevin equation with a classical stochastic 
force accounting for the effect of quantum fluctuations in the field, which causes small fluctuations 
on the particle trajectory. This equation will be useful for the study of stochastic motion of charges 
under the influence of both quantum and classical noise sources, derived either self-consistently (as 
done here) or put in by hand (with warnings). We show the possibility of secular effects from such 
stochastic influences on the trajectory that may impact on the present calculations of gravitational 
waveform templates. 
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I. INTRODUCTION 

Interest in the problem of radiation reaction from par- 
ticle motion in a curved spacetime has seen a rapid in- 
crease in recent years. The back-reaction of emitted radi- 
ation on the particle, known as the self-force, changes the 
particle trajectories (e.g., near a black hole) from a simple 
geodesic motion. The determination of the self force is 
essential to precision calculations of particle trajectories 
and the determination of waveforms from prospective as- 
trophysical sources. Electromagnetic radiation reaction 
in a curved spacetime was first studied by DeWitt and 
Brehm Q. The gravitational radiation reaction equa- 
tion was first obtained by Mino, Sasaki and Tanaka Q 
and Quinn and Wald [j] and others, notably Detweiler 
and Whiting The equation of motion governing a 
scalar charge with radiation reaction was first obtained 
by Quinn Q. For an excellent review, see @- 

Parallel to this there has been detailed work devoted 
to particles and detectors (e.g., atoms) moving in a 
quantum field and in design studies of possible detec- 
tion of Unruh radiation (see, e.g.,0, H, El an d references 
therein). The introduction of worldline path integral 
methods (see and references therein) enable one to 
obtain equations of motion for the charges and the field 
self-consistently. The introduction of open system con- 
cepts and the influence functional method enables one to 
derive stochastic equations with a noise source derived ab 
initio and in a self consistent manner (See 0, ITTl fl2| for 
accelerating detectors, 0,E3>0 for moving charges.) 

In this paper, we use the worldline influence functional 
method to study a particle with a scalar charge moving in 
its own quantum scalar field in a curved spacetime. This 
is a generalization of results obtained in [lOj to curved 
spacetime. We are interested in the radiation it emits 
and its backreaction (radiation reaction) on the trajec- 
tory of the particle and derive the equations of motion for 
the quantum average (expectation value) of the particle's 
position. For those particle trajectory histories which 



become sufficiently decohered, the expectation value be- 
haves classically. The scalar ALD equation we derive in 
this limit checks with the result of Quinn 1]. We then 
include in our consideration the effect of fluctuations in 
the quantum field. We show how it behaves like a clas- 
sical stochastic force and derive a scalar ALD-Langevin 
equation for the particle dynamics with a stochastic com- 
ponent, thus capturing the induced small fluctuations on 
the particle trajectory. This equation will be useful for 
the study of stochastic motion of charges under the influ- 
ence of both quantum or classical noise sources, derived 
either self-consistently (as done here) or put in by hand 
(with warnings). For astrophysical sources with some 
stochastic component this effect may need to be included 
in more accurate calculations of waveform templates. 

In Sec. 2 we describe the worldline influence func- 
tional method and how to obtain the semiclassical and 
stochastic particle dynamics. In Sec. 3 we discuss how 
to regularize the causal Green function in the spirit of ef- 
fective field theory using a quasi-local expansion. In Sec. 
4 we derive the scalar ALD equation. In Sec. 5 we derive 
the stochastic scalar ALD or the scalar ALD-Langevin 
equation. We discuss the non-Markovian nature of the 
noise induced effects and the possibility of secular effects 
from such stochastic influences on the trajectory. In Sec. 
6 we discuss an array of issues pertinent to the present 
problem and approach. In Sec. 7 we summarize our find- 
ings and Appendices A-D provide further details in the 
derivation of certain results given in the text. 



II. RELATIVISTIC PARTICLE-FIELD 
DYNAMICS 



The open quantum system paradigm starts by consid- 
ering a system, or universe, which is partitioned (accord- 
ing to natural physical arguments or some large scale sep- 
aration) into two smaller systems. One subsystem, called 
the system, is assumed to be the one of interest and the 



2 



other, called the environment, contains many more de- 
grees of freedom. An accurate description of the behav- 
ior of the system variables requires knowing the influence 
from the environment due to their mutual interactions. 
A less than accurate description of the overall influence 
of the environment can be obtained by introducing some 
coarse-graining over the environmental variables. 

A. Worldline Influence Functional 

In a coordinate system, assume at some initial time 
ti the quantum statistical state of the combined sys- 
tem S (particle in position zi) and environment E 
(quantum field </?,;) is described by a density matrix 
p(zi, tpi; z[, ifi', ti). In practice, specifying such a state is 
non-trivial since one requires a time-like Killing vector 
to define positive frequency modes and hence a Hilbcrt 
space of states. If the spacetime admits an asymptot- 
ically flat region or is conformally flat then the initial 
state can be constructed. But a general spacetime may 
not admit a time-like Killing vector. Regardless, we will 
sidestep these issues by working at a formal level. At 
some tf > ti the density matrix is evolved to 

= J dz.d^i J dz' i d<p' i K(zf,^f,tf,z i ,(p i ,t i ) 

x pizufc; zi, U) K*(zf,<f>' f , t f ; z[, t f ) 

(2.1) 

where K is the amplitude of the time-evolution opera- 
tor U(tf,ti) = exp{— | J / dt Hs+e[z, 4>]} for the system 
plus environment and has a path integral representation 
given by 

rzf,<t>f 

K(z f ,<f> f ,t f ;z i ,(f> i ,t i )= / VzVc^e^^ (2.2) 
J zi,4>i 

The action describing the system plus environment can 
be written as the actions for the system Ss[z] and envi- 
ronment Se[0] along with an interaction action Si n t[z, <j>] 
between them 

S s+B [z,<l>] = Ss[z} + S E [4>] + S int [z,<l>] (2.3) 

with 



S s [z] = -mo j dr yf-g^u^u" 

Se[4>] = \j d A x^{g^d^d v cj ) ~^ R Rcj ) 2 ) 

S M = Jd* Xj{x ;z M x) 

= dT^HT^{z(T)) (2.4) 



where the current density j(x; z] is given by 

j(x;z}=—^ J dr yf-g^u" 5 ( 2 -5) 

At this stage, r is just a parameter of the worldline and 
not necessarily the proper time. We use an overdot to 
denote differentiation with respect to the worldline pa- 
rameter t. For example, the 4-velocity of a particle in 
Minkowski space (g^ v = rj^ v ) is u a = — z a . Be- 
cause the particle is moving in a gravitational field, one 
should replace ordinary derivatives d/dr by covariant 
derivatives D/dr and so, for instance, the 4-acceleration 
is = u^ p u a and not ii a = = u^d u a . 

To facilitate easier computation it is customary to 
choose the initial density matrix to correspond to a fac- 
torized state of the system and environment. Physically, 
this means that all of the field modes have been uncorre- 
cted with the particle by an instantaneous measurement 
at time ti . Aside from issues about performing this mea- 
surement simultaneously in the spacetime, this choice is 
somewhat unphysical because, as explained in |15l |. an 
infinite amount of energy is required to uncorrelate all 
the modes of the environment (field) from the system 
(particle) at a particular instant of time. For instance, 
in models of quantum Brownian motion with an infinite 
number of environment oscillators the factorized initial 
state results in large transients of the diffusion coeffi- 
cients appearin g in the master equation for the reduced 
density matrix |l§J. The transients appear as a result of 
the high frequency modes of the environment beginning 
to interact and correlate with the system just after the 
initial time. This recorrelation time lasts on the order 
of the inverse of the cut-off frequency used to regulate 
the divergences coming from oscillators of very high fre- 
quencies. For times much longer than this transient time 
the behavior due to the initial factorization is usually 
discounted. Other methods, including the preparation 
function method |l6j| . allow for a somewhat more physi- 
cal initial state by including certain system-environment 
correlations (e.g. system in thermal equilibrium with the 
environment at the initial time), but still seem to suffer 
from some of the problems associated with the factorized 
state 

Having said this we assume that there is a Cauchy 
hypersurface at the initial time ti such that the initial 
density matrix takes the factorized form 

p(z i ,(l)i;z' i ,(p' i ;ti) = ps(*i, z'i'ti) ® PE&i, <P'i,ti) (2-6) 

This simple form eases the manipulations for obtaining 
a description of the reduced particle dynamics. 

After tracing out (a form of coarse-graining) the field 
variables from the density matrix the reduced density ma- 
trix for the system is given by 
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Pr(zf,z' f ;t f ) = / d(j)f p(z f ,<f) f ;z' f7 <f) f ;t f ) 



= / dzidz[ / Vz j Vz' p s {z il z' i \t i )e^ 



|(SsM-S s [*']) 



= dzidz'i Vz Vz 1 ps{z l ,z l l ;t l )eTi { - Ss ^- Ss ^ F[z,z'] 



(2.7) 



The last line introduces the influence functional F[z, z'\, which is given by 

F[z,z'\ = [ d^fdfadti [ 1 V(f> / 0/ I?0 , p £ (^,^;i J )e^ (SB[01+S "' t[ ^ 0] - SEM - S " lt[2 '' ' ]) = e* s w[*>*'] 



(2. 



and Sj„/ is the influence action. In operator language F 
is 

F[z,z'} = Tr B U E+int (t f ,ti;z] p E (U) U E+int (tf,ti;z'] 

(2.9) 

where UE+int (tf,U;z] is the evolution operator that 
evolves the environment variables through its interaction 
with the system. The influence functional can be inter- 
preted as the overlap of the environment states evolved 
forward and backward in time while interacting with dif- 
ferent particle trajectories, or simply as the ensemble av- 
erage of the evolution operator UE+int[z] evolved back- 
ward in time under a different external source z', as can 
be seen by writing l|2.9[) as 

F[Z,Z'] = {U E+mt [z'\ U E+lnt [z]) ens 

= ^2^2pE, aa '(ti) (a'\U E+int [z'} U E+mt [z]\a) 

a a' 

(2.10) 

In the interaction picture, the time evolution operator is 

U E +mt(tf, U; z] = Te^^'^W where the • denotes space- 
time integration so that for two functions, possibly ten- 
sors, A(x a ) and B(x@) 



dx° / d 3 xy^A(x) B(x) (2.11) 



A ■ B 



Assuming that the initial state of the field is Gaussian, 
the influence functional can be calculated exactly giving 



F[z, z'] 



for the influence functional and 



Si E [z, z'\ d = —i h In F[z, z'\ 

= \r- (^ 2 g h ) ■ r + r ■ {i^ 2 G ret ) ■ 3 + 

(2.13) 

for the influence action. The difference and semi-sum 
current densities are defined as 



J = 



J[z\~J[z\ 
j[z] +j[z' 



(2.14) 
(2.15) 



(2.12) 



and the Hadamard Gr and retarded Green's functions 
Gret are 

G H (x, x 1 ) = ({Mx), Mxf)}) - 2 (Mx))(Mx')) 

G ret (x,x') = i0 + (x,£) ([lpl{x),4> I {x')}) 

(2.16) 

The (• • • ) = Tr E p E (- ■ ■ ) denotes quantum expectation 
values in the Gaussian initial state p E of the environment. 
The step function 0+(x, S) appearing in G re t equals one 
in the future of the point x' and zero otherwise. Here, 
£ is a space- like hypersurface containing x' . Had the 
initial state contained non-Gaussian contributions, one 
would have many additional terms involving cubic and 
higher powers of the coupling. Likewise for nonlinear in- 
teractions (e.g. ~ ej[z] -<p n ). Assuming that the coupling- 
is small and that non-Gaussianities are also small one can 
still use (|2.12|l and l|2.13|l as a lowest order approxima- 
tion. 
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B. Semi-Classical Particle Dynamics 

The reduced density matrix for the particle is now writ- 
ten as 

p r (z f ,z' f ;tf) 

= fdz i dz' i J 'vz Jvz' p s ( Zl ,z[-t L )^ ScGEA{ ^ A 

(2.17) 

where the coarse-grained effective action (CGEA) is de- 
fined as 

Scgea[z,z'} 

= Ss{z}-Ss{z / } + S IF [ Z , Z '} 

= S s [z] - S s [z']+r ■ (16^ 2 G ret ) • j+ 



■j~ ■(WTr 2 G H )-f 



(2.18) 



At this point it is worth mentioning that the magnitude 
of the influence functional decays rapidly for two largely 
separated histories since it is Gaussian (to lowest order) 



\F[z,z'}\ = e-TK j ~ 

D 2 



exp 



4fi 
0{zt) 



dr 



dr' ~° SJ 



sr 



Sz a 



(2.19) 

Here, and in the following, z a = z a (r) and z a = z a (r') 
so that an unprimed (primed) index refers to that com- 
ponent of a tensor field or coordinate evaluated at proper 
time t (t') and 



(2.20) 
(2.21) 



The norm of F is equal to the norm of the decoherence 
functional (see below), which is a measure of how much 
the particle's worldline is decohered. So, if the quan- 
tum fluctuations of the field (environment) provide a 
strong enough mechanism for decoherence (this should 
be checked on a case by case basis and will be assumed 
true in the cases under study here) then we are justified 
in expanding the CGEA about the classical trajectory 
z M . Doing so gives 



Scgea[z,z'] = J dr' z°L 



I a ' SScGEA 



5z a ' 



+ 0(z 2 _) 



(2.22) 



Using this expression in the reduced density matrix and 
doing a stationary phase approximation gives the equa- 
tions of motion for the classical worldline z 



SScgea 
Szf(r) 







(2.23) 



Evaluating the functional derivative using, for some test 
function f(x), 



d i x y /^gj (x;z]f(x) 



e ( Du 



47r \ dr 

where wff — g" + u^u 1 ' and u a u c 
gauge), gives 



*+<[z]V,*)/(*) 



m 



dr 



ew^[z\4> ret {z) 



(2.24) 
T (proper-time 

(2.25) 



where 6 r et is the retarded field 



4>ret{x) = e I dr' G ret (x,z a ) 



(2.26) 



The right side of 12.25(1 is the self-force on the parti- 
cle arising from the radiation reaction. Furthermore, 
the vector operator contains the particle's acceler- 
ation implying that the particle moves with an effective 
time-dependent mass equal to m — e(j) ret (z). (See 0] 
for a discussion of evaporating scalar charges based on 
this time-dependent effective mass.) Provided there is 
strong enough decoherence to suppress the quantum fluc- 
tuations of the worldline, the equations of motion for 
the quantum expectation value {z^) are the same as in 
(|2.25l) at tree-level in both the particle and the field ^3 • 
Unfortunatley, (|2.25|) is problematic since the self-force 
diverges on the worldline of the point particle. This di- 
vergence can be traced to the singular coincidence limit 
of the retarded Green's function lim r '_> r G re t — > oo and 
results from the point particle assumption. We will dis- 
cuss how to treat this problem in Sec. IIIII by adopting 
an effective field theory point of view. 



C. Stochastic Semi-Classical Particle Dynamics 

We have made the assumption that the quantum fluc- 
tuations of the worldline are strongly suppressed by the 
decoherence due to the quantum field. Even under strong 
decoherence when the classical trajectory is well-defined, 
the quantum fluctuations of the field can still influence 
the classical motion of the particle through the particle- 
field coupling Si n t- They may show up as classical 
stochastic forces on the particle. In this section we shall 
show how this comes about using the influence functional 
for the reduced density matrix. 

We start by invoking the relation 

(2.27) 
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where N is a normalization factor that is independent 
of the worldline coordinates and £(x) is some auxiliary 
field. Using this relation, the influence functional (|2.12(l 
can be written as 

F[z, z'] = N jv^x) e -it-G5 1 -tHr<-4*t+ie*»Gr«-i+) 

(2.28) 

where appears as an auxiliary function, which can 
be interpreted as a classical stochastic or noise field 0, 
Hsj with an associated Gaussian probability distribution 
functional 

P t [£(x)] = e-te- * 1 * (2.29) 



The fact that this is Gaussian is a direct consequence of 
taking an initial Gaussian state for the quantum field. 
With respect to P%[£\ this implies that £ has zero- mean 
and its correlator is proportional to the Hadamard func- 
tion encoding the information about the fluctuations in 
the quantum field • 

(tix))^ = (2.30) 
({£(*), £(z')}> { - hG H (x,x') (2.31) 

where (...) ( =Nj X>£P ? (...). 

Now the reduced density matrix Ij2.7|l becomes 



p r (z f ,z' f ;t f )=N J dzidzij S Vz JJvz 1 psiz^z'^U) J V^P^] e iS SEA [z,z'; 6 ] 

where the stochastic effective action (SEA) is defined as 

S SEA [z, z'- £] = S s [z] - S s [z'} + j- ■ (-4tt£ + 16n 2 G ret ■ j+) = Sg GEA [z, A - (An) £ ■ f 



(2.32) 



(2.33) 



and Sqqea ^ s the re& l P ar ^ 01 the CGEA. As before, if the worldline is strongly decohered by the quantum fluctuations 
of the environment then we are justified in expanding the SEA around the classical solution 



S S ea[z, z- €i=-j dr'zir],,, \^\ + \ J dr ' J dr " 



, „ ft 2 <?3 



CGEA 



SzrdzT 



■0(z 3 _) 



(2.34) 



Z — Z —Z 



r 



where we have used (|2.23H . expanded j~ 
and defined the stochastic force 



e w. 



[*]£(*) 



in powers of Z- 



(2.35) 



Putting (|2.34(l into the reduced density matrix and using 
the stationary phase approximation gives the stochastic 
equations of motion for the worldline fluctuations z^_ — 



, A2 O-R 

dr'~z a CGEA 



5z^8z a ' 



(2.36) 



It should be emphasized that this equation describes 
the evolution of small perturbations z around the semi- 
classical trajectory that arise from the stochastic man- 
ifestation r] of the quantum field fluctuations. We can 
obtain a stochastic version of the ALD equation (|2.25l) if 
we add the equations of motion for the classical dynamics 



0D CGEA 







to the left side. Then we may write 



0D CGEA 



6zt 



(2.37) 



(2.38) 



Evaluating this functional derivative gives the stochas- 
tic semi-classical particle dynamics for the full worldline 

Z[l % [l i z^ 

Du 

m = eWu[z]<j> ret {z) + n^[z] (2.39) 

Notice that both the deterministic and the stochastic 
components of the self-force can push the particle away 
from its geodesic motion with respect to a fixed back- 
ground spacetime. However, we should keep in mind that 
this equation is really only valid to linear order in the 
fluctuations z since higher order terms will correspond 
to quantum corrections that we have been neglecting. Of 
course, we could have obtained l|2.39|l by using a station- 
ary phase approximation in the reduced density matrix 
(|2.32|) without having expanded around the classical tra- 
jectory. However, we believe that the validity of (|2.39() 
to linear order in the fluctuations would not have been 
as transparent. 

The stochastic correlation functions of the force ij^ can 
be evaluated using the knowledge of the £ correlators 
above. As was commented in the previous paragraph we 
must evaluate these correlation functions along the clas- 
sical trajectory z to be consistent with the linearization. 
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The mean of the stochastic force is zero 

(»7 M [*]> f = e^[z] (t(z)) ( = (2.40) 

and the symmetric two-point function of the stochastic 
force is 

({nn[z]>iln'[z]})e = he2 W{n [z]w,j,i)[z] G H (z a ,z a ) 

(2.41) 

This shows that the noise r/^ is multiplicative, colored 
and depends on the particle's initial conditions through 
the classical trajectory. The noise correlator also generi- 
cally depends on the field's initial conditions as is seen by 
the appearance of Gh in the equation. The Hadamard 
function Gh does not vanish on a space-like hypersurface 
implying that the quantum correlations in the environ- 
ment are non-local. However, the noise correlator above 
is evaluated on time-like separated points only. For equal 
proper times r' = r the Hadamard function diverges so a 
suitable regularization procedure must be used in order 
to make sense of (I2.41f) near coincidence. 

In the next two sections we describe a new method in- 
troduced in for regulating the ultraviolet divergences 
in the retarded Green's function. We introduce a high 
energy scale in the quantum field below which the low 
energy point particle dynamics is expected to be valid. 
For long times as compared to the inverse of the high en- 
ergy scale, we invoke a quasi-local expansion in order to 
obtain the relevant contributions to the particle's motion. 
Our work is a generalization of 10] to curved spacetime. 
We follow Poisson's description of scalar radiation reac- 
tion in curved spacetime and use units where c = G = 1 . 
(We change the particle physics notation in [lj} (+, -, -, 
-) to the MTW notation (-, +, +, +) of |fj.) 

III. REGULARIZATION OF THE RETARDED 
GREEN'S FUNCTION 

Before we discuss regularizing the retarded Green's 
function it will be beneficial to introduce some notation. 
Following Q we define the step function 9 + (x, E) to be 1 
for any point x to the future of a space-like hypersurface 
E and otherwise. This is a generalization of the ordi- 
nary step function to curved spacetime. Using this, one 
can define the following distributions 

S+(a(x,x')) = 9 + (x,V)6(a(x,x')) (3.1) 
9+(-a(x,x')) = 9 + (x,T,)9(-a(x 1 x')) (3.2) 

where a(x,x') is Synge's world function along the 
(unique) geodesic linking x = z(t) and x' = z(t'). Tak- 
ing the lightcone centered on x', with x' on the space- like 
hypersurface E, it follows that 5+(c) has support along 
the forward lightcone only while 9 + (— a) is one in the 
causal future of x' and vanishes everywhere else. 

In (|2.1(il) the retarded Green's function is calculated 
from the commutator of the interaction picture field </>/. 



Recall that <fii evolves under the free dynamics so that 

(□-6^)07 = (3.3) 

and so is the homogeneous solution to the full field equa- 
tion in the Heisenberg picture (satisfying the same initial 
conditions) 

{n-£ R R)<t> H =j[z} (3.4) 

where <f>n = 4>i + 1 G re t ■ j [z] for a classical source. Be- 
cause the source is classical only the first term in the so- 
lution gives a contribution to the commutator and hence 
to the retarded Green's function 

Gret ~ ([$h(x),$ H (x')]) = ( [fa(x), (3.5) 

Finally, since the commutator is state-independent it fol- 
lows that G re t for the quantum field 4>h and for the cor- 
responding classical field <p are the same. 

A. Hadamard Expansion 

This allows us to invoke Hadamard's ansatz for which 
the retarded Green's function for a scalar field in a curved 
3 + 1 dimensional spacetime has the form Q 

Gret {X, X ) 

= A 1 / 2 (x,x')6+(a(x,x')) + V(x,x') 9 + (-a(x,x')) 

(3.6) 

which is the sum of a "direct" part (proportional to 5+) 
and a "tail" part (proportional to + (—<r)). Notice that 
the appearance of 5 + and 9 + ensures that G re t has the 
correct causal structure. In order for this ansatz to be 
valid, the points x and x' must be connected by a unique 
geodesic. Otherwise, ambiguities arise from the appear- 
ance of caustics when parallel propagating a tensor field 
from x 1 to x or vice versa. This will be assumed in the 
following. The function A(x, x') is the van Vleck determi- 
nant and the function V(x, x') satisfies the homogeneous 
Klein-Gordon equation in curved spacetime 

(n g -UR)V(x,x') = Q (3.7) 

with boundary data determined by the restriction of V 
to the forward lightcone, denoted by V. The restriction 
of V is found by solving 

v> Q + \ « - 2) v = \ {(□, - UR) Al/2 }_ 

(3.8) 

Knowing all of the light-like geodesies emanating into 
the future from x', one can integrate this equation to 
construct V on the light-cone and hence V inside the 
light-cone. Of course, this just solves for V(x, x') for any 
x to the causal future of x' and must be solved again for 
different values of x' . This makes determining V very 
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difficult and tedious for generic spacetimes. However, if 
the spacetime possesses enough symmetry (e.g. de Sit- 
ter, some FRW cosmologies) then V can be constructed 
relatively easily 0] . 

From the mean and stochastic equations of motion 
()2.25(l and (|2.39|1 it is necessary to compute the restric- 
tion of G re t to the particle's worldlinc 

G ret (z a , z a ') = A 1 / V, z a ') 5+(a(z a , z a ')) 

+ V(z a ,z a ')9+(-a(z a ,z a ')) (3.9) 

Using the Hadamard ansatz and l|2.26[l it is easy to see 
that for two (time-like separated) points on the particle 
worldline, 5+(cr) has support only at coincidence, when 
t' = r. The contribution from the direct term is then 

J ' dr' A^ 2 {z a , z a ') 9+(z a ,S) 6(a(z a ,z a ')) 
= £ dr' A 1 / 2 (z a ,z a ')S(T' -t) (^j 1 (3.10) 

This expression is ultraviolet divergent since da/dr eval- 
uates to zero for t' — r and so the retarded field is not 
well-defined (in fact, it is infinite) when evaluated at the 
location of the point particle. This feature is not limited 
to just a quantum field as it appears even for a classical 
field since they share the same retarded Green's function. 

The presence of this divergence suggests that a reg- 
ularization procedure is needed to render the equations 
of motion l|2.25|l and l|2.39|l finite. Several procedures 
have been proposed in the literature for doing this. In 
[H, , the retarded field is evaluated near the worldline 
so that the divergence can be tracked in the limit that 
one approaches the worldline. This divergence renormal- 
izes the particle's mass. Writing the retarded field in 
terms of functions singular (f>g and regular 4>n on the 
worldlinc allows, after calculating their gradient and ex- 
panding them near the worldline, for the self-force to 
be evaluated unambiguously along the worldlinc. This 
method can also be used for electromagnetic and gravi- 
tational self-force calculations. In another method devel- 
oped in |20L l2ll l22l I2.I |2 1| , a mode-sum regularization 
procedure is used in which the self-force is decomposed 
in terms of the multipole moments of the field. From 
this one subtracts the moments of the singular part of 
the self-force and resums over the multipoles to obtain 
a self-force that is finite on the worldline. In [25| the 
"massive field approach" is used in which an auxiliary 
massive field is introduced to subtract away the ultravio- 
let divergence from the massless field. This procedure is 
similar to Pauli-Villars type regularization and has also 
been used in |l3j ] to regularize the direct term in G re t ■ 

Other regularization methods have been used for the 
radiation reaction of non-relativistic particles without 
gravity present j2||. Notable of these include the 
extended-charge models of [23 in which the particle is 



described as a rigid volume with a non-zero size. Un- 
fortunately, this type of prescription cannot be extended 
to relativistic models since rigid bodies are incompati- 
ble with relativity ^j. Furthermore, the successes of 
quantum field theory (e.g. as applied to the scattering of 
point particles), has proved to be an adequate descrip- 
tion of the low-energy effective dynamics of point par- 
ticles to ~TeV scales. It seems unnecessary, then, to 
introduce non-trivial structure to the particle to study 
the self-force. 



B. Regularization motivated by Effective Field 
Theory 

With this consideration in mind we introduce yet an- 
other regularization for the self- force motivated by [To| . 
All theories for the fundamental interactions (e.g. QED, 
electroweak theory, QCD) can be interpreted as effective 
theories, in that they are are low energy limits of a more 
complete theory valid at higher energies, yet they pro- 
vide an excellent description of low energy phenomena, 
at least for renormalizable theories, without requiring the 
details of the complete theory. For energies much higher 
than the energy scale of the effective theory, new physics 
is likely to become important. 

In this spirit we introduce a regulator A for the field. 
Above the energy scale of this regulator, E\, new physics 
is assumed to occur. For energies much lower than E\, 
the dynamics (|2.25() and i|2.39[l will be sufficiently accu- 
rately described by using a regulated quantum field. This 
approach has been taken in [lOj in deriving the Abraham- 
Lorenz-Dirac (ALD) equations in flat spacetime. It is 
easily extended to motions in a curved spacetime since 
the direct part of the retarded Green's function is a local 
quantity on the worldline. 

Regularization is achieved here by choosing any suit- 
ably smooth function to approximate 6+ for large values 
of the regulator A. Following [l(| we replace 5+ with 

S+(a(z a ,z a ')) -» e+(-a(z a ,z a '))J^A 2 e- 2Ai<r2 ^ 

= e + (-a(z a ,z a '))gi t (z a ,z a ') (3.11) 

In the limit of infinitely large A we have that 8+ g^ et ap- 
proaches 6+ . The integral of the direct part of G ret (I3.1U|) 
receives a contribution from <5 + at coincidence implying 
that a = 0. The function g^ et is smooth but approxi- 
mates S + well only if A 2 a 3> 1. This approximation will 
not hold if a is identically zero. Nevertheless, we will 
assume that a is small and approaching zero but A is 
such that A 2 a 3> 1 still holds. This separation of scales 
allows us to do a quasi-local expansion below in which 
the self-force will be expanded near coincidence. 

Making the replacement l|3.11|l in G re ± gives the regu- 
lated Green's function 
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Gt t {z a ,z a ') = 9 + (-*(z a 1 z a ')){A 1 ^(z a ,z a ')gUz a ,z a ') + V(z a ,z a ')}=9 + (-a){A 1 ^gi t + v} (3.12) 

where an abbreviated notation has been used on the far right side to ease the appearance of later expressions. Also, 
it should be clear that in the limit A — > oo that G^ et — > G re ± ■ 
Putting G A ei into (£Z5) gives 

mo^=e 2 V*«) ^ T/ dr^(r-rO{AV^A t + F } (3 .13) 

where the quantities in the integral are evaluated along the classical trajectory z. Passing the derivative through the 
integral gives 

[Vp (t - r')] + jT dr> V 5 , (AV^A t + J 

where the [. . .] denotes the coincidence limit r 1 — > r of the quantity it contains. This notation looks the same when 
denoting a functional of a function (e.g. w^z]) but the context should be clear. From Appendix A one can show 
that 

[V^(r - t')\ = - [V„(a a u a )} = -u v (3.14) 

and since iu„"[z] projects vectors onto a direction orthogonal to the mean 4- velocity then [V^ (r — r')] gives no 

contribution leaving 




(3.15) 



= e 2 ^ 4>ret{i) + e 2 w»[z] jT dr> { (v 2 . A'/ 2 ) g^ t + A 1 / 2 V 2 , ffr A et + V^v} (3.16) 

I 



to ■ 



dr 



At 



4>ret{z) + e 2 W"[z] 



A 1/2 5 A et 



r 



It will be convenient for later manipulations to define 
the tail term as the integral over the past history of the 
gradient of V, 

cf^ il (z(T))=ej T dr'V^V (3.17) 

The tail term will turn out to be responsible for the non- 
Markovian and history-dependent dynamics of the parti- 
cle's evolution. Consequently, this presents a significant 
source of difficulty when trying to determine the parti- 
cle's motion since one has to know the entire past history 
of the particle's worldline in order to make predictions 
about its current position and speed. 

IV. EFFECT OF RADIATION REACTION ON 
PARTICLE TRAJECTORY 

In order to simplify some of the difficulties of solving 
the nonlinear integro-differential equation (|3.1(i|) we uti- 



lize the effective theory viewpoint discussed in the previ- 
ous section. Since <? A et is strongly peaked around a = 
then the contributions to the self-force will be largest 
when t' ~ t. We therefore introduce a quasi-local ex- 
pansion by expanding 6 + A 1 / 2 g^ et around coincidence 
t' = t. This yields a well-defined description of the self- 
force in terms of relevant and irrelevant quantities. The 
relevant terms are those that either renormalize certain 
parameters of the particle or otherwise affect the parti- 
cle dynamics when A goes to infinity while the irrelevant 
terms are those that give no contribution in the same 
limit. 



A. Quasilocal Expansion 

We begin by expanding the van Vleck determinant and 
<? A et for small values of a(z a , z a ). Following ||| the ex- 
pansion of the square root of the van Vleck determinant 
around er, a a = is 



A 1 / 2 (z a ,z a ')^l + -R^(z) a a a? 



1 

36 



R a0 , 7 a a o p a 1 



(4.1) 
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where the . . . denotes terms of higher order in the expansion. The covariant derivative with respect to the worldlinc 
coordinate of the above quantity is 

V 2 , A 1 / 2 ^, z a ') = 1 R a0 (z) a% + ^ Rcp-A*) <r a ^ + ... (4-2) 
Lastly, cr"^ can be expanded around small a to give 

<J a . v = g a v -\R a lu5 {z)^a 5 + (4.3) 

so that the gradient of A 1 / 2 becomes 

V 2 , A 1 / 2 (f° , Z"' ) = ~ i?„ Q ( Z >°-i R^af3u;\ (z) ^ ^ + . . . (4.4) 

Next, we expand these quantities and g^ et around r 1 = r + s with s such that As ^> 1. The expansions of some 
relevant quantities involving the world function a are (see Appendix A for details) 

s = r' -t (4.5) 

a(r,T') = - l -s 2 + 0(s i ) (4.6) 

/ x s 2 DuJt) s 3 D 2 u u (t) ... , , , 

Using these results, the function g^ et is 

= ^ A 2 C 2 A 4 cr 2 = /J A 2 e -AV/J + 0(s 6 } (4 g) 

Lastly, we expand Vf<7^ e4 m powers of s 

y-^9ret ~ ^9ret - ) ~ \Mt) + - + - +0( S )f— (4.9) 

Applying these results to (|3.16|) gives 

e / \ 9 /x o I _™'x \ iy a 



(m - y 3(1) (r) + e 2 c (0) (r) - e 2 jf rir' U(z Q , z Q ') 



fir 

= ew;[z] C a {z) + e 2 J2 (s ( „+i)(r>i" +1) M + c (b) (»>W [z]) (4.10) 
n=l 

with r — t — Tj being the elapsed proper time since Tj. 

The expressions for the r-dependent coefficients (<?(«) and ci n y) and the trajectory-dependent vectors (uffl and 
can be found in Appendix B. However, only the relevant factors in the limit of very large A will be given here 



A fl r 4 A 4 \ , , . 2 7 /4 / 5 r 4 A 4 

«(0)(r) = -^174^7 (j, — j ff(i)(r) = -A 7 ^, 



2 



(4.11) 



1 H ^ 4 A 4 \ / x 2/3 r 4 A 4 \ n . 

c « w = a 7 b —J 9(2)(r) = t^U* — J (4 - 12) 

= ^;[f]iC«a (4-13) 
n^[z\ = i S„ u&>[z] = I ^ (4.14) 
I 

where 7(0, b) = T(a) — L(a, b) is the incomplete gamma scale of ~ A -1 after which they are effectively constant, 
function. The coefficients (I4.11|) - (|4.12|) vary over a time- 
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C (0)\ 





rA 



FIG. 1: Time-dependence of the coefficients appearing in 
14. Ill - 14. 121 . The functions C(o) and have been divided 
through by A so that they can be displayed along with 
and 3( 3 ). 



mass mo gets shifted by an infinite amount 5m = 
e 2 ((7(!)/2 — C( )). The mass is then renormalized to 
m ren = mo — 5m, which is a constant with respect to 
r in the infinite A limit. If the regulator is finite but very 
large then the shift 8m is finite and large and the dressed 
mass has a very weak dependence on time (see Fig. Q). 

If we now make A infinite and use the limiting values 
for 5(2) and c (1) 



Urn g(2) (r) = 1 = lim c (1) (r) 



(4.17) 



for a constant elapsed proper-time r then the low-energy 
effective particle dynamics is 



m re „(r; z] 



m, 



dT ~» V ' J " 

where the renormalized effective mass is 



F^(r) + U[z] (4.18) 



W i Z\ — Tfl T 



e 2 / dT'V(z a ,z a ') (4.19) 



Figure 0) shows the r-dependence of these functions. 
Note also that, at the initial time, all of these functions 
vanish. More will be said below concerning these prop- 
erties and their implication for the validity of the quasi- 
local expansion. 



B. Scalar ALD Equation 

Most of the terms on the right side of l|4.10|l are irrele- 
vant in the sense that they are inversely proportional to 
powers of A so that for time-scales much longer than A -1 
these terms can be ignored. In this limit the effect of the 
high energy physics is ignorable and the effective particle 
dynamics is described by 

( -i Du » 

+ J c (i)M <K Rua(z) u a + ew»[z] <Pl ail (z) 
+ C(A- 1 ) (4.15) 

where the time-dependent and trajectory-dependent ef- 
fective mass is 

m(r; z] 

= m -e 2 Q 5( i)(r) - c (0) (r) +^ dr'V(z a ,z a ') 

(4.16) 

and an external force F^ xt responsible for accelerating 
the charge has been included. 

In the limit that A — > oo, the terms involving C( )(r) 
and 5(i) (r) both diverge linearly with A and the bare 



and the self-force is 

= 2 Da 



e -^+ e -R„ a {z)u a + e^\z) 



(4.20) 



Equations (|4.18|l through i|4.2U|l are the main results of 
this section. The appearance of the tail term (|3.17|) im- 
plies that the radiation reaction is nonlocal because of its 
dependence on the past behavior of the classical world- 
line. This equation for the radiation reaction from scalar 
charges was first obtained by Quinn based on earlier 
work of Quinn and Wald Q and Mino, Sasaki and Tanaka 
on gravitational radiation reaction. 

Detweiler and Whiting [j| have obtained this result by 
decomposing the derivative of the retarded field </> r et into 
a singular piece </>^, containing the diverging contribu- 
tion that renormalizes the mass, and a regular piece <j>^, 
which contributes to the self-force and is regular on the 
worldline. From the previous section we can construct 
these quantities within our regularization scheme to find 
that for A — > oo 



"' : ^ A2s " + e (^ (1)W " C(0)(r) ) ^dV 
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R{z)u lJ 



(4.21) 



Notice that the first term of the singular part does not 
contribute to renormalizing any physical parameters (at 
the level of the equations of motion) since </>^ is projected 
onto a direction orthogonal to u^. Likewise, the first term 
of the regular part does not contribute to the self-force 
on the particle. 

In a flat spacetime, the tail term vanishes since there is 
no curvature available to focus the radiation emitted in 
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the past onto the particle at the present. Furthermore, 
(|4.18|) reduces to the ALD equation for a scalar field 

e 2 

m ren ^ = F™\t) + j w^[z] Q„ (4.22) 

This equation was derived in the open quantum system 
formalism in fijjl. 



V. LINEARIZED FLUCTUATIONS AROUND 
THE CLASSICAL PARTICLE TRAJECTORY 

We now study the effects of the quantum field fluctu- 
ations (as classical stochastic forces) on the low energy 
dynamics of the particle. Instead of working directly with 



the equations for the linearized fluctuations (|2.36() , which 
contain the singular functional derivatives, we begin with 
(|2.39l) and assume that the retarded field has been regu- 
larized by the short-distance regulator A in a quasi-local 
expansion. 

A. Stochastic ALD Equation 

For large finite values of A i|2.39[) becomes the ALD- 
Langevin equation 

Dn 

™<r,z\-^ = K 'fW + UW + iti-nz] (5.i) 

where the (regulated) self-force is 



fM 



Slz] (j9(2)(r) ^ + ^c (1) (r)fl m (z) U a +e rw)+0(A ( - 1) ) (5.2) 



Of course, one must remember that these expressions are only valid up to linear order in the fluctuations z about 
the mean worldline z since higher orders correspond to quantum corrections that we have assumed to be negligible. 
Expanding l|5.1|l in orders of the fluctuations, the time-dependent mass and self-force are given by 

m(r;z] = m(r;z]-e 2 ^ dr' ' ~z v ' jp J dr" V(z a , z a " ) + 0{S 2 ) (5.3) 

fM = fM + J dr' z-'-^fM"} + 0(S 2 ) (5.4) 
Calculating the functional derivative in the mass equation gives 

m(r;z] = m(r;z] - - e 2 f dr' z" 'V ? y V \z a , z a ') + O '(5 2 ) (5.5) 

where the factor of 1/2 in the second term comes from evaluating a delta function on the upper boundary of the r' 
integral. Notice that the linear terms in z vanish in flat space-time since V is identically zero for a massless scalar 
field in 3+1 dimensions. 

Simplifying the self-force fluctuations is slightly more involved. The calculation amounts to performing the varia- 
tional derivative on /„ but keeping in mind to expand out the covariant derivatives, which depend on the classical 
worldline coordinates. The result is 

fM] = fM] - K MA ~ za - iMA z a - ™»M\ * a + r MA ' z ' a + o(AM (5.6) 

The time- and trajectory-dependent coefficients are complicated expressions that are not particularly illuminating. 
These are recorded in Appendix C. 

Combining the linearized mass and self-force into l|5.1|l and using the fact that z satisfies the mean equations of 
motion l|2.23|l results in 

M m [z] S a + IVJz] z a + KMA z a - e2 ^l dT ' z u V^V(z a , z a ') = r m [z] 2° + vM + 0(AM (5.7) 



where and 
K^M] = K lia [z]-m(T-,z]dg a T^uptr' --^^'(2) 
Miia[z] = ff MQ m(T;z] +to mq [z] (5.8) 



W z l = iM z ] - 2m ( T ; z ] r ^( Q u /3) ( 5 - 9 ) 
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The dynamical equation (|5.7(l for the fluctuations about 
the classical particle trajectory is the main result of this 
section. This is a linear integro-differential equation for z 
with a third derivative term and contains time-dependent 
coefficients that depend on the non-Markovian behavior 
of the mean trajectory. Furthermore, because of the inte- 
gration over past times the last term on the left side de- 
pends on the history of the fluctuations. Notice that this 
term vanishes in a flat background so that the fluctua- 
tions then obey a third-order differential equation, which 
is Markovian in the sense that given a mean trajectory 
the fluctuations do not depend on their own past history. 
So, given a solution to the mean equation of motion for 
the classical worldline one could, in principle, solve for 
the fluctuations induced by the quantum field fluctua- 
tions. 

The notation in (|5.7|l has been chosen suggestively 
since the left side resembles a damped simple harmonic 
oscillator with time-dependent mass, damping factor, 
and time- and history-dependent spring constant. No- 
tice also that the "effective mass" M^ a [z] is not diag- 
onal implying that the inertia of the fluctuations be- 
haves differently in different directions. This feature is 
exhibited in all of the other coefficients (T^q,, K^ a and 
r^a) and suggests that the fluctuations of the trajectory 
in one direction are linked with the fluctuations in the 
other space-time directions. Also, from the expression 
for the radiation reaction on the fluctuations r^ a [z] it 
should be noted that this is explicitly independent of the 
background curvature and effectively projects z „ onto a 
direction orthogonal to the mean 4- velocity u v . 

If the stochastic term r/^ is ignored then (|5.7|) describes 
the evolution of small perturbations away from the semi- 
classical trajectory and so is useful for studying the linear 
response of the trajectory to small perturbations away 
from the mean worldline. In other words, setting r/^ to 
zero (|5.7|) gives the linearization of the ALD-Langevin 
equation around the semiclassical worldline z. General- 
izing this to the self-force due to linearized metric per- 
turbations could be useful for investigating the stability 
of numerical calculations of the inspiral of a small mass 
black hole into a large mass black hole, for example [30| . 

B. Memory and Secular Effects 

A particularly interesting feature related to this is the 
effect of the non-Markovian term appearing in (|5.7|) 

-e 2 a^J dr' z v 'V- z ^V{z a ,z a ') (5.10) 

If the fluctuations z grow then the linearization may not 
be valid for all times and so one should then include quan- 
tum corrections to our stochastic semi-classical equations 
if the environment is a quantum field as is assumed here. 
The growth of the fluctuations could be significantly in- 
fluenced by the past behavior of V V V and the history 
of the fluctuations. The behavior of this term on the 



fluctuations is not known yet because of the difficulty in 
solving integrodifferential equations and in computing V 
for many spacetimes. However, a study in de Sitter space 
(where V is known |rill9j]) may provide a simple arena for 
understanding the role of the non-Markovian term and 
the nature of the solutions to (|5.7|) . 

Aside from these technical considerations, these equa- 
tions are applicable for any type of noise on the particle 
trajectories. In many cases, the source of the noise act- 
ing on the system of interest may not be known and is 
put in by hand in a phenomenological description of the 
particle dynamics via the stochastic equation 

Dn 

Mr; A-^ = *T M + //.(*) + 4 ( 5 - n ) 

where the superscript (+) is to remind us that this term 
was added in by hand, as opposed to being derived, like 
our treatment of r]^[z\. This stochastic force could have 
a classical origin (e.g. high-temperature thermal fluctua- 
tions of surrounding gas) or it could have no known single 
identifiable origin. Furthermore, since the noise 77+ is not 
derived from an initially closed system it is likely to be 
inconsistent with the dynamics of the trajectory. (See, 
e.g., 0) In any case, one would also have to specify the 
noise correlator (r)£ (t)^, {r')) n + as it suits the model. 

With this proviso (no guarantee for consistency) the 
analysis of this section carries over. Given any kind of 
noise the fluctuations around the mean trajectory of the 
particle moving through its own (classical) field subjected 
to the self-force from radiation reaction is given by (|5.7|l , 
but with r)fj,[z] replaced by 77+ 

M m [z] l a + r MQ [z] ~z a + K^ a [z] ~z a 

-e 2 a J dr' fV f ,-F(z a ,z a ') 

= r, a [zYz' a +r ] ;+0(A- 1 ) (5.12) 

Since in this discussion we don't need to worry about 
quantum corrections from higher order loops in the ef- 
fective action we could go beyond the linear order in the 
fluctuations of the particle trajectory and expand the so- 
lutions to l|5.11|) in powers of the coupling constant e ( 
denoted in the subscript) z = zq + z\ + Z2 + . . .. As- 
suming that 77+ is 0(e) and depends on z and recalling 
that is quadratic in the coupling we find a non-local 
and causal contribution to the total force on the particle 
coming from the fluctuations of the stochastic force (see 
Appendix D for details) 

FfT ift = J dr' dr" F^(r,r',T")( V ;(zS')vi(zf)) v+ 

(5.13) 

which is of the same order as the self-force. It seems 
that the stochastic noise would cause the particle to drift 
off from the background trajectory Z( ) determined by 
the external force. The deviation could build up over 
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time as indicated by the integral above. This may have 
interesting consequences for astrophysical sources with 
some stochastic behavior described by a classical noise. 
If such physical situations exist, this noise- induced drift 
may give rise to a secular effect which could alter the 
waveform templates (of events expected to be seen by 
gravitational interferometers like LIGO and LISA) cal- 
culated without including such stochastic secular effects. 

A similar expression to (|5.13l) can be derived in our 
open quantum system framework by writing the effective 
action to cubic order in the fluctuations. Doing this re- 
veals a term like ()5.13f) but this requires a much more 
careful analysis that goes beyond the stochastic semi- 
classical approximation adopted here. 



VI. DISCUSSIONS 

A. The Quantum Regime and the Validity of the 
Quasi-Local Expansion and Order-Reduction 

The previous sections focused on the Feynman- Vernon 
formalism and various approximations to obtain the low- 
energy effective dynamics of the particle, both for its 
semiclassical (mean) and stochastic semi-classical (mean 
and stochastic fluctuations) motion. Here, the domain 
of validity of the quasi-local expansion and this semi- 
classical treatment will be explored and compared with 
the relevant scales for weak and strong radiation damp- 
ing. 

In the effective field theory paradigm a regulator A is 
introduced for controlling the ultraviolet divergences ap- 
pearing in the direct part of the retarded Green's function 
such that AV > 1 with a small and approaching zero. 
After expanding a near the coincidence limit the time- 
scale of the quasi-local expansion At = s is governed 
by 

At>A~ : (6.1) 

Recall that for times larger than ~ A -1 the time- 
dependent coefficients in (|4.15|l and JO} rapidly ap- 
proach constant values (see Fig|I). But there are other 
scales to consider. We have been working at tree-level in 
both the particle and the field so that At must be much 
larger than the scale for 1-loop field and 1-loop parti- 
cle corrections. This implies that the time interval At 
should be much longer than the time-scale for creating 
pairs 

At > — = A c (6.2) 

TO 

where Ac is the Compton wavelength associated with 
the scalar particle. The dynamics l|4.15|l and l|5.1|l are 
therefore valid using this framework provided that 

At > Ac > A" 1 (6.3) 



Nevertheless, the presence of the third T-derivative in 
(|4.15l) , H5.1fl and (|5.7|) requires the specification of the ini- 
tial position, velocity, and acceleration to obtain unique 
solutions. This is problematic, since for a vanishing ex- 
ternal force F^ xt one still requires an initial acceleration 
to solve the equations. But, if there is no force accel- 
erating the charge then what causes the particle to ac- 
celerate? Furthermore, the particle may experience un- 
bounded acceleration so that its kinetic energy increases 
with time to infinity. Since scalar dynamics is an energy- 
conserving theory where then could this energy arise? 
These problems are related to the infinite self-energy of 
a (classical) charged point particle. 

In the classical theory of scalar fields and particles 
these issues can be resolved if one instead gives the par- 
ticle a finite size ro for which the self-energy is roughly 
e 2 /ro 31]. If this energy composes its rest mass then 

e 2 

r (6.4) 

TO 

and represents the "size" of the particle. In j^l), an 
approximation, amounting to an asymptotic expansion 
in powers of ro, called the Landau approximation (also 
known as order-reduction), is employed to obtain solu- 
tions that require only an initial position and velocity 
and are also free from run- away solutions. The Landau 
approximation converts the ALD equation (of third or- 
der) to the so-called Landau-Lifshitz equation (of second 
order). We will use these names to distinguish between 
these equations. 

In order-reduction, the lowest order solution is found 
by simply ignoring the self-force so that the radiation 
damping is assumed weak. The time-scale of the dynam- 
ics is then determined mostly by the external force so 
that if Ffi xt varies on a scale \ ex t then At ~ X ex t- In 
curved spacetimes the self-force will be weak if ro <C At 
and the length scale associated with the spacetime curva- 
ture Xr is large (i.e. small curvature) so that A# 3> ro. It 
then follows that for weak damping in the semi-classical 
domain that the quasi-local expansion and the Landau 
approximation are valid provided that 

At - \ ext > r , Ac > A" 1 and \ R > r (6.5) 

Recently, in the context of plasma physics, [2(| has 
investigated the validity of the Landau approximation 
for the classical ALD equation by numerically integrat- 
ing the Landau-Lifshitz equation forward in time and, 
using the final position, velocity, and acceleration from 
that, integrating the ALD equation backward in time. 
If the initial position and velocity of the particle differ 
significantly from the backward-evolved solution of the 
ALD equation at the initial time then one can assume 
the Landau approximation has broken down. Koga does 
this for a counter-propagating electron and ultraintcnsc 
laser beam (intensity ~ 10 22 W/cm 2 ). He finds that the 
Landau approximation is valid so long as the laser wave- 
length Aq is greater than the Compton wavelength. For 
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Ao much smaller than Ac, he finds disagreement between 
the solutions of the Landau-Lifshitz and ALD equations. 
However, these equations cannot be fully trusted since 
quantum effects should become important. In this do- 
main, while we cannot directly apply our results to this 
problem we can use the closely related closed-time-path 
(CTP) formalism to incorporate the effects of the quan- 
tum loop corrections to the (quantum) particle dynamics. 

B. Decoherence 

This brings us to an important issue, namely, how 
strongly must the particle worldline decohere in order for 
a stochastic semi-classical treatment to be applicable? A 
measure of the decoherence of a system is the decoher- 
ence functional D[a(z), a(z')] where ot[z) specifies a 
particular coarse-grained history of the worldline. If we 
take the history to be a(z) — z then the decoherence 
functional is 

D[a{z),a{z')] = D[z,z'\ (6.6) 
= p s (zi,4)ei Saa **l'>'1 (6.7) 

and is closely related to the CTP effective action [33| . 
In fact, the norm of D is proportional to the norm of 
the influence functional. The system is said to deco- 
here if \D\ approaches zero. We showed earlier that 
the CGEA is proportional to and so there is mini- 
mal decoherence for two nearby histories. This implies 
that coarse-graining the quantum field must be comple- 
mented by an additional coarse-graining of the worldline 
in order for the particle to be sufficiently decohered. If 
the mechanism for decoherence is efficient then this addi- 
tional coarse-graining should be minimal and the stochas- 
tic semi-classical description that we have used in this 
paper is applicable. For our problem, recall that the 
CGEA contains an imaginiary part that is proportional 
to j~ ■ Gh ■ j~ so that the norm of D is 

\D[z,z']\ = e -TKr-^G H )-r ( 6g ) 

showing that the quantum field fluctuations in the envi- 
ronment is the mechanism for decoherence. It is therefore 
necessary to understand the behavior of the Hadamard 
function to show that there exists a well-defined stochas- 
tic semi-classical regime for the particle. This requires 
knowing the field modes on the background space-time, 
which is a significant problem in its own right. For simple 
scenarios the decoherence functional can be calculated 
and the decoherence time approximated. But, in general, 
one needs to determine if decoherence is fast enough on 
a case by case basis. 

C. Problems with putting noise in by hand 

Towards the end of Sec. 5 we replaced the classical 
stochastic manifestation of the quantum fluctuations by 



a noise source 77+ that has some specified noise correla- 
tor. Aside from the issue about keeping only to tree-level 
in the particle and field variables, how is using 77+ differ- 
ent from all of the sophisticated machinery used in the 
previous sections? 

Unlike in our approach, inserting a source of noise by 
hand implies that it can have any correlator one wants. 
The stochastic two-point function of 77^ 

({vi(r),vt(r')}) v+ =N^(r,r'), (6.9) 

otherwise known as the noise kernel, can be chosen at 
will. Some physical reasoning, for instance, a high tem- 
perature environment, might suggest that N IjIU (t,t') ~ 
S(t — t'), i.e. white noise. While white noise certainly 
simplifies the calculations, there is little justification for 
simply stating the form of the noise kernel and expecting 
the dynamics of the system to be consistent with that 
source. Using the Feynman- Vernon formalism, coarse- 
graining over the quantum environment naturally gives 
rise to the noise kernel appropriate for the stochastic force 
fluctuations 

<{run^r']}> € 

= he 2 w (p [z a ]w u) [z a '\ G H (z a , z a ') (6.10) 

This coarse-graining ensures that the dynamics of the 
system evolves consistently with that of the environment. 
This is an important statement that cannot be overstated 
because it is this feature which gives rise to consistent 
fluctuation-dissipation relations (FDR). Inserting noise 
by hand most likely violates such a relation. 

While FDR's won't be discussed here in great de- 
tail, have studied the FDR's for n particles inter- 
acting with a quantum field for certain trajectories (e.g. 
uniformly accelerated). They have shown the existence 
and self-consistency of these relations and further intro- 
duce correlation-propagation relations relating the corre- 
lations of particles to each other through the (causal) 
propagation of the quantum field. These correlation- 
propagation relations can only be formed if the parti- 
cles interact directly through causal influences. The open 
quantum system viewpoint, therefore, shows the interre- 
lated influences of the system and environment through 
these types of self-consistent relations. This is an impor- 
tant feature lacking in a model that has noise put in by 
hand. 

Another problem with adding noise in by hand is that 
the noise can be a significant contribution to the quantum 
two-point function of the system if one is interested in the 
quantum particle behavior. So, having the correct form 
for the noise kernel is important for getting the correct 
expression for the quantum two-point functions. 

It has been shown in using the CTP generating 
functional for quantum correlation functions of the sys- 
tem variables, for quantum Brownian motion as well as 
for stochastic semi-classical gravity (see next section), 
that certain quantum two-point functions are related to 
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stochastic two-point functions. As an example, the sym- 
metrized quantum two-point function of a Brownian par- 
ticle's position can be written as 

i ({£(*), £(*')}) = ((X{t)X(t')) i ) Xipi (6.11) 

where X is a solution of the QBM Langevin equation 
L ■ X = £ for the appropriate linear operator L(t, t'), and 

I 



where TV = (^(^1)^(^2))? is the noise kernel, given by the 
stochastic correlator of £. The first term involving the 
homogeneous solution of the Langevin equation repre- 
sents the dispersion in the initial conditions and is called 
the intrinsic fluctuations. The second term involving the 
noise kernel, and hence the quantum fluctuations of the 
environment, represents the correlations with the envi- 
ronment through the stochastic fluctuations and is called 
the induced fluctuations. It turns out that if the homo- 
geneous solution Xq decays exponentially fast then for 
late times, at least times larger than the decay time, the 
quantum two-point function is determined entirely by the 
noise kernel, that is, the induced fluctuations. And so, in 
this sense, all of the information about the quantum cor- 
relations of the system degrees of freedom is encoded in 
the stochastic correlations. This shows another drawback 
to putting noise in by hand. The induced fluctuations 
contain information about the quantum fluctuations of 
the system variables. But this will not be true with some 
arbitrary noise kernel chosen at will. 



D. Similarities with Stochastic Semiclassical 
Gravity 

The features of the particle dynamics seen in the above 
discussions are typical of nonequilibrium open quantum 
systems. History-dependent behavior is present in the 
equations of motion for the system and if a renormaliza- 
tion procedure is required it is usually a time-dependent 
prescription, as seen earlier with the renormalized mass 
"Ven (■?")■ Furthermore, the noise correlator Uti.lUfl is 
generically non-local in time and is determined by the 
quantum fluctuations of the environment variables. This 
formalism does not allow for arbitrary noise kernels since 
this would destroy the self-consistency between the sys- 



(• ■ -)xi,pi denotes the average over all possible initial po- 
sitions and momenta with respect to the reduced Wigner 
function for the initial state of the particle. The solutions 
of the QBM Langevin equation consist of a homogeneous 
part, containing all the information about the initial con- 
ditions of the system, and a term describing effects due 
to the interactions between the particle and oscillators. 



(6.12) 



(6.13) 



tern and environment evolution. A particular example 
that contains these features is stochastic semiclassical 
gravity, which we will briefly describe and compare with 
below. 

Stochastic semiclassical gravity (SSG) is a self- 
consistent theory of the stochastic dynamics of a clas- 
sical spacetime containing quantum matter fields. SSG 
goes beyond semiclassical gravity, for which the geom- 
etry is driven by the expectation of the (renormalized) 
stress tensor, in that the quantum field fluctuations also 
contribute to the spacetime dynamics through a classical 
stochastic source. The spacetime is therefore driven by 
both the quantum expectation value of the renormalized 
stress tensor and a classical stochastic stress-tensor-like 
object, For an introduction and review of this sub- 
ject see [34| and [35j for a discussion of the domain of 
validity of SSG. 

As an open quantum system, the quantum field fluc- 
tuations are coarse-grained using the CTP formalism of 
Schwinger and Keldysh (SK) to study the self-consistent 
evolution of the (classical) geometry. The quantum fluc- 
tuations manifest themselves as stochastic noise thereby 
imparting a stochastic nature to the spacetime. The 
resulting Einstein-Langevin equation for the linearized 
metric perturubations h a b is 

Gab \9 + h] = n( T« [g + h] + k U[ 9 ] (6.14) 

The superscript W denotes that those quantities con- 
tain all terms to first order in the metric fluctuations 
h a b. It should be noted that the finite parts of the coun- 
terterms needed to cancel the divergences coming from 
the stress tensor expectation value have been absorbed 
into the definition of (T^) ren . The renormalized stress 
tensor expectation value (evaluated in a Gaussian state) 
contains an integration over the past history of the met- 
ric fluctuations and so the dynamics is generally non- 



X(t)=X (t)+ / dt' G ret (t,t')Z(t') 



The symmetrized quantum two-point function becomes 

±({x(t),x(t')}) = (X (t)X (t')) Xupi + I dh I A 2 G re4 (Mi)A^i,i 2 )G ret (t 2 ,0 



r 
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Markovian. This is like what is seen in the ALD-Langevin 
equation (|5.1|l where the tail term <fi t J L ' ll (z) is analogous 
to the expectation value of the renormalized stress tensor 
in (|6.14|) . The (covariantly conserved) stochastic source 
tensor £ ab has zero mean and its correlator is given in 
terms of the Hadamard function of the stress tensor fluc- 
tuations t ab = f ab - (T ab ) 

({Lb(x;g],£,cd(x';g}}) i = ft ({i ab {x; g],i cd (x' ; g]}) 

(6.15) 

The correlator of the stress tensor fluctuations on the 
right side does not vanish on a spatial hypersurface. This 
reflects the fact that the quantum field correlations are 
themselves non-local. Compare this with the correlator 
in 12.31fl which is also non-local. 

SSG also suffers from runaway solutions since the fi- 
nite contributions to the counterterms needed to cancel 
the divergences appearing from the expectation value of 
the stress tensor are quadratic in the curvature. This 
makes SSG a theory with derivatives higher than two, 
similar to the ALD type equations derived above, which 
were of third order in the r derivatives. A version of the 
Landau approximation can be used to reduce the order 
of the Einstein-Langevin equation to two thereby yield- 
ing well-behaved solutions free of the pathologies typical 
of higher-order derivative theories. Of course, one needs 
to be careful to use the Landau approximation at scales 
that are consistent with the derivation of the Einstein- 
Langevin equation. 

Finally, the symmetrized quantum two-point functions 
of the metric fluctuations h ab can be written in terms of 
intrinsic fluctuations, representing the dispersion in the 
initial conditions, and induced fluctuations, encoding the 
information about the fluctuations of the quantum mat- 
ter |35l ]. Just like with the particle motion, one can- 
not simply use any noise kernel for modeling stochastic 
metric fluctuations. One needs to do a careful analysis 
that ensures the self-consistency of the metric and quan- 
tum matter dynamics and the existence of fluctuation- 
dissipation relations. 



VII. SUMMARY 

In this paper we have derived the scalar ALD equa- 
tion for the quantum expectation value of the world- 
line for a scalar charged point particle interacting with 
its own quantum field as it moves in a curved space- 
time. Our equation l|4.18|l for the low-energy effective 
particle dynamics agrees with the results obtained ear- 
lier by 0, U @, 0| . If the quantum fluctuations in the 
field strongly decohere the worldline then we can ignore 
the particle's quantum fluctuations and obtain the semi- 
classical motion (|4.18|l . Invoking an effective field theory 
point of view, the singular behavior of the field's retarded 
Green's function can be regulated. For sufficiently short 
times At = s but still long enough At A -1 compared 



with the inverse cutoff frequency A a quasi-local expan- 
sion can be used to obtain the contributions relevant to 
the self-force and those that are irrelevant in the infinite 
A limit. This renormalizes the mass of the particle and 
shows explicitly the appearance of the expected Da^/dr 
term that is characteristic of radiation reaction. 

While the time-dependent coefficients l|4.11ll - (|4.12ll ap- 
pearing in (|4.15l) seem to suggest that only the initial po- 
sition and velocity of the particle are needed, and hence 
that a resolution of the problems of pre-acceleration and 
run-away solutions has been reached, one should keep in 
mind that the quasi-local approximation breaks down for 
short time intervals, which includes the instant at the ini- 
tial time. A more careful analysis would need to include 
a more physical initial state than the factorized one used 
here. 

Fluctuations in the quantum field is expected to affect 
the particle's motion causing it to fluctuate by an amount 
z around the mean trajectory z given by solutions to 
the semi-classical equation (|4.15|l . We derived such a 
stochastic force ^^[z] with correlators from the fluctua- 
tions of the quantum field and a scalar ALD-Langevin 
equation 1)5.10 . The dynamics of z 1)5. 7|) contains a non- 
Markovian contribution through the past history of the 
particle fluctuations. Depending on the behavior of this 
term the fluctuations might grow to be large indicating 
a breakdown of the linear approximation. In that case it 
requires the inclusion of quantum corrections in order to 
follow the nonlinear evolution of z. On the other hand, 
ignoring the noise altogether in (|5.7|l . one can test the 
stability of numerical simulations in inspiral studies, for 
example |3Cj . 

Instead of the noise rj^ derived here from fluctuations 
of quantum fields one can replace it with some other clas- 
sical noise 77+ suitably chosen to model some stochastic 
source in a phenomenological description. We can still 
use (|5.7|) to study the effect of such noises on the parti- 
cle trajectory fluctuations. However, since the origin for 
the noise is no longer due to a quantum field we need 
not worry about keeping up to linear order in z in the 
ALD-Langevin equation. Instead, expanding the solution 
for small coupling constant e and taking the stochas- 
tic expectation value shows that there is, in general, a 
non-vanishing force (|5 . 131) coming completely from the 
correlations of the stochastic force. Along with the self- 
force, this noise-induced term would cause the particle to 
drift off of its background trajectory determined by the 
external force (or off of its geodesic motion if F^ xt van- 
ishes). We hope to explore the consequences of this noise- 
induced drift in an astrophysical setting and find observ- 
able effects on the waveforms of the radiation emitted by 
the particle and detected by gravitational interferometers 
like LIGO and LISA. 

In Paper II we will apply the same techniques here 
to study the self-forces and the stochastic semi-classical 
motions of electric charges coupled to an electromagnetic 
field and of small black holes coupled to the background 
spacetime of a massive black hole, respectively. We hope 
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to apply these results to more physical situations such 
as the motion of charges in strong external fields and 
gravitational radiation reaction. 
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APPENDIX A: GEODESIC COORDINATES AND 
THE QUASI-LOCAL EXPANSION 



FIG. 2: For two points on a worldline 7 one can construct a 
unique time-like geodesic j3 connecting them. 



In this Appendix we derive the quasi-local expansion 
of a a in 1|4.7|) and some of the relations appearing in 
Section IV for computing the quasi-local expansion of the 
regulated direct part of the retarded Green's function. 

We begin by invoking the so-called geod esic coordi- 
nates along and near the worldline [3ll l36| . These co- 
ordinates will be denoted with a hat. In geodesic coor- 
dinates the metric tensor and the connection coefficients 
are constructed to vanish along the entirety of the world- 
line 







pA I 







(Al) 



so that a vector field on the worldline behaves like a vec- 
tor field in flat spacetime. The displacement from A to 
B in Fig© can be written as 



A ~ a de f ~c 
Az = z 



Using a Taylor series to express z a in terms of z c 
its derivatives results in 



(A2) 
and 



Az a = 



2! dr 
s 2 Du a 
2! 



s a u 
6! dr 2 ~ 
s 3 D 2 u a 
' 61 



dr 6! dr 2 

where, in the last line, we have used the fact that the 
components of the connection vanish in this coordinate 
system. 

Flatness along the worldline implies that the tangent 
spaces at each point on 7 can be identified with each 
other so that tensor manipulations along 7 can be done 
in a single tangent space. This implies that Az a is a vec- 
tor in the tangent space that we can interpret as a dis- 
placement vector at B. Transforming back to the original 
coordinates gives 



Az a = 



dx a 
OS? 

-su° 



AZ? 
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s 3 Da a 
6\ "dr" 



(A4) 



Using geodesic coordinates appropriate for the 
geodesic connecting A and B (denoted with a bar instead 
of a hat) one finds that 



a = (A - A )t = y -y 



(A5) 



which represents the coordinate difference between A and 
B in the tangent space at B. Here we take i a to be the 
tangent vector to the geodesic (3 connecting A and B with 
affine parameter A" and with coordinates y such that 
V\a = and y\s = In the original coordinates 







(A6) 



Since Az a and a a both represent the (coordinate) differ- 
ence between A and B in the same tangent space then 
they are related by the coordinate transformation from 
the hat to the bar coordinates 



f^^AP 



(A7) 



(A3) This immediately implies that 



a a (z a ,z a ) 



a s 2 a s 3 Da a 

-su -a 

2! 6! dr 

^^.d^ U M 

n— 1 



(A8) 



There is an interesting interpretation one may attach to 
this result. Pick a point x' in the normal convex neigh- 
borhood of x in some spacetime. These points are con- 
nected by a unique geodesic described by a a . If these 
two points intersect any worldline describing the motion 
of some particle then in terms of the particle's velocity 
and its derivatives a a — Az a . This equality therefore 
describes a (many-to-one) mapping between particle mo- 
tions and geodesies. 
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From the identity 2a = a a a a it is possible to show 
that 



9(n)(r) = - 



„2 4 

2 24 V ; 



(A9) 



2 (n+6)/4 / „ r 4 A 

liW A 7 ( 1+ 4'T 



This also provides a relation between the parametriza- 

tion of the geodesic and the proper time of the worldline ,. , , 

° r Ine corresponding worldlme-depcndent vectors 

through 



arc 1 



(A-A') 2 = s 2 + ^+0( S 5 ) 



(A10) 



since 2a = — (A — A') 2 is the amount of elapsed proper 
time along the geodesic. 

Using these results, the covariant derivative of a scalar 
function, say g^ et (a) , is 



da\ 1 dg^ et (a) 



ds 



^% A £i+0(s 4 ) (An) 
S OS 



APPENDIX B: COEFFICIENTS AND VECTORS 
IN REGULATED DYNAMICS 

The r-dependent coefficients (r = t — r.j) appearing 
in the equation for the regulated semi-classical dynamics 
are 



C (n)(r) 



ds 

o n 



9(™- 1 )/ 4 f\ 4- n r 4 A 4 



J 



(0) 


= 






(1) 


1 

2 






(2) 


1 

3 


< 


_Det Q 
dr 


(3) 


1 

4 


< 


D 2 a 
dr 2 



■ a fJ ,R al 3U a u 13 



v (0) 


= -a M 




1 

= — w 
6 




1 

= — w 
6 



-R 7 u a vP 



(B2) 

(B3) 
(B4) 

(B5) 

(B6) 

(B7) 
(B8) 

(B9) 



For finite and large A the first-order correction to the 
self- force f a [z] is 



(g m {r)uf\z\+c {2) {r)v^[z_ 



(BIO) 



(Bl) which in flat spacetime reduces to ^ 9<3)(t) w£ [z] D d ^ !) 



APPENDIX C: COEFFICIENTS IN LINEARIZED STOCHASTIC DYNAMICS 

The time- and mean worldline-dependent coefficients appearing in l|5.7J) as a result of linearizing the self-force 
are 



K m {T]z] = yS(J)(r) ^ ( 

"< [2] (y C(i) (r)u" V 2 « ^ + eV z -« # aiZ (*)) 



7i ? 



(CI) 



p 2 /„2 rj= 2 

--^[f]c (1) (r)i2, Q - 2tt<"0 / , )a f - ff(2) (r)^ + - c (1) (r)i^ + e<fr U {z) 



(C2) 



2 

( T ; f ] = y 5(2) (r) u^" [z] (^2 r^u^ g 7)Q + vPVJ v g^ 



(C3) 
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r m (r;z] = —g(2){r)w lJ 



(C4) 



APPENDIX D: NOISE-INDUCED DRIFT 

In this section we outline the details for obtaining the 
noise-induced force on the particle l|5.13|l . We begin with 
the ALD-Langevin equation with noise added in by hand 



dr 



(Dl) 



If we expand the worldline coordinates in powers of the 
coupling constant e 



(D2) 



where the subscript denotes the order of the expansion, 
then to lowest order one finds 



(D3) 



where D^/dr is the covariant r-derivative defined with 
respect to Zq . If the external force vanishes then this is 
nothing more than the geodesic equation. 

The first order correction to the trajectory can be 
found by solving 



dr 



which has the formal solution 



*?(t)= / dr'gZ^r')^^') 



(D4) 



(D5) 



and g^ et is the retarded Green's function for z\. 
To second order in e we find that z 2 satisfies 



m ren {^^- ~ 2rf jU (2 )wo Q u 2/3) - zldqV^{z Q )u^uo^j 

dr' V(z%, z£') + /„[*„] + z^v+izo) 



rpext 



en J n 



+ mren \T ^ {z )u^u w + 2zjd z ^ V ^ (z )u[ a u fj) + ^zlz^d z yd z sT^(z )u^u ^j 



(D6) 

Taking the stochastic average we find that the terms after the self-force give a non-vanishing contribution to the force 
on the particle 



Ft ift (r) = ^ 1 V z ,r 1 +(z )+m ren [T^mp + 2zJd zZ T^u o0) + ^Jzld^T^u^ 

Substituting in for zi we find that F*"*/* can be written in terms of the noise kernel associated with 77+ 

Ft ft (r) = J 'dr' C'V 2S «( Z „ A )e(4')), + 

+ m rm J dr' dr" j ^ d^d^T^z^Uop g^ t g S r a e 't + 2 ^ r f M ( 2 o ) 9ret u a(/3 ~f- 9rlt 



(D7) 



Notice that the first term in (|D7(I comes from linearizing the stochastic force whereas the other terms have their 
origins in the kinetic term m[z\a^. Finally, we arrive at (|5.13() by defining the kernel Fff as the integrand appearing 
above 



FfT^r) = J dr'dr"F^(r,r',r")«(4')^(^o")>, + 

(D8) 
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